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LONG-TERM  GOALS 

The  goal  is  to  investigate,  theoretically  and  through  analyzing  existing  data,  the  role  that  sea  spray 
plays  in  transferring  heat,  moisture,  and  momentum  across  the  air-sea  interface,  especially  in  high 
winds.  Ultimately,  I  plan  to  develop  simple,  fast,  physics-based  parameterizations  for  these  air-sea 
fluxes  for  use  in  large-scale  models,  especially  those  simulating  tropical  and  extratropical  storms. 

OBJECTIVES 

1 .  By  analyzing  existing  data  sets  and  CBLAST  flux  data  as  they  become  available  from  the 
hurricane  experiments,  develop  simple,  fast  parameterizations  for  the  air-sea  sensible  and  latent  heat 
fluxes,  the  total  enthalpy  flux,  and  the  surface  stress  in  high  winds,  where  sea  spray  is  mediating  all  of 
these  exchanges. 

2.  Theoretically  extend  these  parameterizations  to  high  winds,  up  to  hurricane  strength  (—60  m/s). 

3 .  Collaborate  with  large-scale  modelers  to  implement  and  test  these  fonnulations  in  state-of-the- 
art  coupled  atmosphere-ocean  models. 

APPROACH 

This  work  is  theoretical  and  analytical;  it  has  no  experimental  component.  Andreas  is  the  only  CRREL 
participant,  aside  from  an  occasional  undergraduate  student;  but  he  has  been  collaborating  with  large- 
scale  modelers  elsewhere — primarily  Will  Perrie  at  Bedford  Institute,  Dartmouth,  Nova  Scotia — to 
implement  his  spray  parameterization  in  mesoscale  stonn  simulations. 

Microphysical  theory  establishes  how  rapidly  spray  droplets  can  exchange  heat  and  moisture  in  a  given 
environment.  Theory  also  predicts  how  sea  spray  production  should  depend  on  wind  speed.  The 
analytical  part  involves  developing  parameterizations  for  the  various  spray  transfer  processes  by 
simplifying  model  results  or  by  synthesizing  various  data  sets  and  observations.  Checking  the 
parameterizations  against  available  data  is  also  another  aspect  of  what  I  call  analytical  work. 

Conceptually,  momentum  and  sensible  and  latent  heat  can  cross  the  air-sea  interface  by  two  routes:  as 
interfacial  fluxes  that  are  adequately  parameterized  by  the  COARE  bulk  flux  algorithm  (Fairall  et  al. 
1996)  and  as  fluxes  mediated  by  sea  spray.  In  low  winds,  say  10  m/s  or  less,  the  spray  route  is 
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negligible.  As  spray  concentration  increases  with  increasing  wind  speed,  however,  the  spray  route 
becomes  increasingly  important  until,  in  hurricane-strength  winds,  it  is  the  dominant  air-sea  exchange 
route 

(e.g.,  Emanuel  2003;  Andreas  2004a).  I  have  been  developing  a  unified  algorithm  that  predicts  the 
flux  contributions  through  both  routes. 


My  algorithm  assumes  that  the  total  sensible  heat  (Hs  Tot),  latent  heat  (HL_Tol),  and  momentum  (tto1) 
fluxes  that  would  be  measured  by  eddy-correlation  instruments  at  a  height  above  the  spray  layer  are 
the  sums  of  the  interfacial  [i.e.,  Hs,  HL,  Tmt(0)]  and  spray  [i.e.,  Qs.sp,  Ql.sp,  tsp(0)]  contributions: 


^s.Tot  Hs  Qs.sp  ’ 

(la) 

HL,Tot  =  HL  +  Ql.sp  ’ 

(lb) 

Tot  =  Tjnt  (0)  +  Tsp  (0)  . 

(lc) 

In  (lc),  the  zeros  remind  us  that  the  total  momentum  flux  imposed  by  the  large-scale  atmospheric 
pressure  gradient  comprises  both  interfacial  and  spray  components  right  at  the  sea  surface.  In  a 
modeling  sense,  the  total  fluxes  represented  as  the  left  sides  of  (1)  would  serve  as  the  lower  flux 
boundary  condition  for  an  atmospheric  model. 


My  flux  algorithm  uses  the  COARE  version  2.6  algorithm  (Fairall  et  al.  1996;  Andreas  and  DeCosmo 
2002)  to  compute  Hs,  HL,  and  Tint(0)  in  (1).  I  estimate  xsp(0)  in  (lc)  from  a  simple  mechanistic  model 
and  knowledge  of  the  rate  at  which  sea  spray  is  generated  (Andreas  2004a).  Finally,  a  full 
microphysical  model  that  I  have  developed  (Andreas  1989)  can  estimate  Qs,sp  and  Ql,sp,  but  this  model 
is  too  computer-intensive  for  large-scale  modeling.  Hence,  I  have  greatly  parameterized  the  results 
from  that  model  such  that  my  current  flux  algorithm  predicts  Qs,sp  and  Ql,sp  from 


Q 


S,sp 


Pw  cw  (Ts  Teq,100 )  Vs  (u, )  , 


Ql,SP  =  PwLvj1 


r('Cf,5o) 

50  pm 


r  VL  (u«)  • 


(2a) 

(2b) 


In  these,  pw  is  the  density  of  seawater;  cw,  the  specific  heat  of  seawater;  Lv,  the  latent  heat  of 
vaporization  of  water;  Ts,  the  sea  surface  temperature;  and  Teq,ioo,  the  equilibrium  temperature  of  spray 
droplets  that  originally  formed  with  a  radius  of  100  pm. 
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Furthermore,  in  (2b),  r  is  the  radius,  as  a  function  of  time  (t),  of  a  droplet  that  started  with  a  radius  of 
50  pm.  In  general,  for  any  spray  droplet, 

r(0  =  req  +  (r0  -  req)eXP(-t/Tr)  »  (3) 

where  ro  is  the  initial  radius  of  the  droplet,  req  is  its  equilibrium  radius,  and  xr  is  the  e-folding  time  for 
the  droplet’s  evolution  to  this  equilibrium  radius.  In  (2b),  Xf^o  is  the  atmospheric  residence  time  of  a 
droplet  that  started  with  a  radius  of  50  pm.  Thus,  in  essence,  I  approximate  Ql,sp  from  the  behavior  of 
50-pm  droplets.  Likewise,  (2a)  estimates  Qs,sp  from  the  behavior  of  100-pm  droplets. 

Finally,  in  (2),  Vs(u„)  and  VL(u„)  are  empirical  functions  of  the  friction  velocity,  u*  =  (xTot  /  pa  )’/2 , 
where  pa  is  the  air  density.  I  showed  these  function  in  my  last  annual  report  (also  Andreas  2004b). 

WORK  COMPLETED 

Equations  (2)  and  (3)  contain  three  microphysical  quantities:  Teq,  req,  and  xr.  For  completeness,  I  also 
must  mention  that  a  spray  droplet’s  temperature  evolution  can  be  characterized  by  a  fourth 
micrometeorological  quantity,  xT,  the  e-folding  time  for  the  droplet  to  cool  from  Ts  to  Teq  (e.g., 

Andreas  1990).  That  is,  a  very  good  approximation  for  a  spray  droplet’s  temperature  evolution  is 

=  exp(-t/xT)  .  (4) 

As  Aeq 

My  full  microphysical  model  (Andreas  1989)  computes  all  four  of  these  quantities  by  numerically 
tracking  the  evolution  of  a  spray  droplet.  But  again,  such  tedious  computing  is  too  time-consuming  to 
make  (2)  viable  in  large-scale  modeling.  Hence,  during  this  past  year,  I  developed  approximation 
formulas  and  wrote  FORTRAN  code  to  compute  Teq,  req,  xT,  and  xr  for  a  wide  range  of  droplet  sizes 
and  marine  conditions  (Andreas  2005). 

Equation  (2b)  includes  another  time  scale,  Xf,  a  spray  droplet’s  atmospheric  residence  time.  I  estimate 
this  as  (Andreas  1992) 


xf 


h1/3 

2uf  (r„)  ’ 


(5) 


where  H1/3  is  the  significant  wave  height,  and  Uf(ro)  is  the  terminal  fall  speed  of  a  droplet  with  initial 
radius  ro.  [Hence,  in  (2.6),  xr.50  is  evaluated  from  (5)  for  ro  =  50  pm.]  The  relation  I  have  been  using  to 
predict  H1/3,  however,  seems  to  give  unusually  high  waves;  Xf  is,  therefore,  too  long;  and  the  latent  heat 
flux  computed  from  (3)  and  (2b)  is,  in  turn,  too  big.  This  year,  with  the  help  of  a  Dartmouth 
undergraduate  student,  I  therefore  completed  the  analysis  of  18  years  of  wave  data  from  NOAA  ocean 
buoys  and,  thereby,  developed  a  new  parameterization  for  H1/3  (Andreas  and  Wang  2006). 
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Lastly,  this  year,  Perrie  et  al.  (2005)  published  simulations  of  several  midlatitude  Atlantic  storms  that 
incorporated  my  spray  flux  algorithm.  Our  sensitivity  studies  suggest,  in  general,  that  accounting  for 
spray  heat  transfer  produces  better  simulations  of  maximum  storm  winds  and  minimum  central 
pressure  than  simulations  that  include  only  the  usual  interfacial  heat  fluxes. 

RESULTS 
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Figure  1.  Temperature  and  radius  evolution  of  a 
sea  spray  droplet  (initial  salinity  34  psu)  that 
started  with  a  radius  of  100  pm  and  a 
temperature  of  20°C.  The  droplet  is  ejected  into 
air  with  a  temperature  of  18°C  and  a  relative 
humidity  of  90%.  The  dashed  lines  show  the 
results  from  my  full  microphysical  model;  the 
solid  lines  represent  the  exponential 
approximations  (3)  and  (4),  which  also  rely  on 
my  approximation  formulas  for  Teq,  tt,  req,  and 
z,.  The  temperature  curves  are  basically 
inseparable  out  to  20  seconds,  while  the 
exponential  approximation  predicts  a  slightly 
slower  radius  evolution  than  the  full  model. 
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Figure  2.  Significant  wave  height  versus  10-m 
wind  speed  from  a  year  (2003)  of  hourly  data 
from  NOAA  buoy  41001,  which  was  in  4427  m  of 
water  150  nautical  miles  east  of  Cape  Hatter  as. 
The  yellow  symbols  are  the  hourly  data.  The 
black  circles  are  averages  for  wind  speed  bins 
that  are  0.5  m/s  wide;  the  error  bars  are  ±2 
standard  deviations  in  the  bin  means.  The 
curves  show  our  new  model,  (7),  and  the 
traditional  wave  height  relation,  (6).  The 
traditional  prediction  for  H1/3  clearly  gives  waves 
that  are  much  too  high  for  wind  speeds 
above  10  m/s. 


The  approximation  formulas  that  I  developed  for  predicting  Teq,  req,  tx,  and  xr  work  well.  With  these, 
Teq  can  be  quickly  estimated  to  within  0.02°C,  req  and  xT  can  be  estimated  to  within  about  5%,  and  xr 
can  be  estimated  to  within  about  10%  for  typical  oceanic  conditions.  Figure  1  supports  this  statement 
by  showing  computations  for  the  evolution  of  a  1 00- pm  droplet  based  on  my  full  microphysical  model 
and  on  the  exponential  approximations  (3)  and  (4),  with  my  new  formulas  for  finding  Teq,  req,  xt, 
and  xr. 
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The  results  of  our  analysis  of  the  wave  data  from  the  NOAA  ocean  buoys  were  surprising.  My  spray 
flux  algorithm  currently  uses  a  typical  equilibrium  wave  model  to  estimate  H1/3  in  (5)  (e.g.,  Earle  1979; 
Andreas  1992;  Taylor  and  Yelland  2001): 


H1/3  =  0.030  U/0  ,  (6) 

which  gives  H1/3  in  meters  for  the  wind  measured  at  10  m,  U10,  in  m/s.  Figure  2,  however,  shows  that 
this  relation  does  poorly  in  representing  a  year  of  hourly  data  from  one  of  the  buoys  we  studied. 

Figure  2  is  typical  of  18  years  of  data  that  we  analyzed  from  12  different  buoys  off  the  northeast  coast 
of  the  U.S.  in  waters  with  depths  ranging  from  18.9  to  4427  m.  Each  buoy  record  exhibits  a  range 
where  the  average  wave  height  is  constant  for  10-m  winds  up  to  4  m/s.  Above  this  range,  H1/3 
increases  as  Uf0  ,  but  not  nearly  as  rapidly  as  (6)  suggests. 

From  this  analysis,  we  have  developed  a  scheme  to  predict  H1/3  for  all  water  depths;  it  takes  the 
general  form 


H1/3  =  C(D)  for  0<UIO<4m/s,  (7a) 

H1/3  =  a(D)  +  b(D)Uf0  for  4m/s<U10  .  (7b) 

Here,  H1/3  is  again  in  meters  when  U10  is  in  m/s;  and  C,  a,  and  b  are  empirical  functions  of  the  water 
depth  (D)  only.  The  curve  labeled  “New  Model”  in  Figure  2  is  based  on  (7).  For  the  18  years  of  data 
that  we  have  analyzed,  (7)  predicts  the  measured  hourly  wave  height  at  each  buoy  with  a  correlation 
coefficient  that  is  typically  0.7  or  better. 

IMPACT/APPLICATIONS 

The  unified  turbulent  flux  algorithm  that  I  have  developed — represented  formally  as  ( 1 )— (3 ) — has 
three  features  that  are  not  all  present  in  any  other  air-sea  flux  algorithm:  It  explicitly  recognizes  two 
routes  by  which  heat  and  momentum  cross  the  air-sea  interface,  the  usual  interfacial  route  and  the 
spray-mediated  route;  it  has  been  verified  against  data  (Andreas  2003,  2004b);  and  it  is  theoretically 
based  and,  therefore,  can  be  extrapolated  to  high-wind  conditions,  where  we  currently  need  an  air-sea 
flux  algorithm  but  have  few  reliable  data  on  which  to  base  an  empirical  one. 

Developing  such  a  high-wind  air-sea  flux  algorithm  is  one  of  the  central  objectives  of  CBFAST.  But 
we  still  need  to  see  if  such  an  algorithm  improves  predictions  of  ocean  stonns.  I  have  been  trying  to 
answer  that  question,  primarily  through  my  collaboration  with  colleagues  at  Bedford  Institute.  Our 
simulations  with  an  earlier  version  of  my  flux  algorithm  (Version  1.1;  Andreas  2003)  suggest  that 
including  the  spray  heat  fluxes  in  a  mesoscale  atmospheric  model  gives  better  predictions  of  the 
intensity  of  extratropical  storms  than  does  a  more  conventional  surface  flux  parameterization  when 
central  pressure  and  maximum  surface-level  wind  speed  are  used  as  metrics  for  stonn  intensity  (Fi  et 
al.  2003;  Perrie  et  al.  2004a,  2004b,  2005). 
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Although  I  developed  approximation  fonnulas  for  the  microphysical  quantities  Teq,  req,  Tt,  and  xr 
specifically  for  use  in  my  air-sea  flux  algorithm,  the  need  for  such  microphysical  quantities  arises  in 
many  other  applications.  For  instance,  marine  aerosols  continually  evolve,  and  these  formulas  can 
quickly  predict  the  rates  of  these  changes.  Although  my  formulas  currently  treat  only  saline  droplets, 
the  microphysical  theory  on  which  they  are  based  applies  to  all  aqueous  solution  droplets. 
Consequently,  with  minor  changes  to  account  for  the  chemical  properties  of  the  solute,  my  formulas 
could  be  adapted  for  cloud  droplets  that  have  formed  on  a  variety  of  cloud  condensation  nuclei — not 
just  on  sea-salt  particles. 

TRANSITIONS 

Besides  the  journal  articles  and  conference  papers  that  I  have  written  to  describe  my  work  on  sea  spray 
and  the  resulting  bulk  flux  algorithm,  I  have  developed  an  electronic  “kit”  that  contains  the  instructions 
and  FORTRAN  tools  necessary  to  implement  version  2.0  of  this  algorithm.  I  have  distributed  this  kit 
to  several  CBLAST  collaborators — namely,  Bill  Frank  at  Penn  State  and  Shouping  Wang  at 
NRL-Monterey — and  to  several  others. 

The  transition  of  my  work  that  has  progressed  the  furthest,  however,  is  at  Bedford  Institute,  where  Will 
Perrie  and  his  colleagues  have  introduced  my  unified  surface  flux  parameterization  into  the  Canadian 
mesoscale  compressible  community  model  (MC2)  and  have  been  doing  simulations  of  Atlantic  storms 
with  it.  This  work  is  already  documented  in  several  papers  (Li  et  al.  2003;  Perrie  et  al.  2004a,  2004b, 
2005). 

RELATED  PROJECTS 

I  have  no  other  support  for  this  type  of  work. 
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